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We study numerically the real-time history of the cosmo- 
logical electroweak phase transition, as it may take place in 
the Standard Model or in MSSM for mn ^ rnw according to 
recent lattice results. We follow the nucleated bubbles from 
the initial stages of acceleration and rapid growth, through 
collisions with compression waves resulting in slowing down 
and reheating to Tc, until the final stages of slow growth and 
evaporation. We find that collisions with compression waves 
may make the bubble walls oscillate in the radial direction, 
and that reheating to Tc takes generically place. 

PACS numbers: 98.80. Cq, 47.75. +f, 95.30.Lz 

The electroweak phase transition in the early Universe 
may have had important consequences, like baryogene- 
sis [Q] . These depend crucially on the non-equilibrium de- 
tails of the first-order transition ||] . The general features 
of the real-time history of first-order cosmological phase 
transitions have been known for quite a long time 
and, e.g., the properties of bubble growth soon after the 
nuclcation period have been studied in detail [4-6]. The 
whole real-time history has also been computed in some 
approximation but a more detailed hydrodynami- 

cal investigation has so far been missing. In addition, 
the relevant parameter values are reliably known only 
after recent lattice simulations The purpose of this 
paper is to present a detailed hydrodynamical study of 
the complete real-time history of the electroweak phase 
transition for realistic parameter values. 

Before the transition the universe is in the high- 
temperature 'symmetric' phase. As the temperature falls 
below the critical temperature Tc, the symmetric phase 
becomes metastable. Bubbles of the 'broken-symmetry' 
phase nucleate at r„ < Tc, and grow rapidly. There are 
two modes of bubble growth at this stage: detonations, 
which are supersonic, and defiagrations, which are usu- 
ally subsonic |l^,|l^. In a deflagration the bubble wall 
is preceded by a compression wave, whereas in a deto- 
nation it is followed by a rarefaction wave. When the 
bubbles collide, the evolution becomes complicated. We 
find that the compression waves from neighboring bub- 
bles may cause large-amplitude back-and- forth motion of 
the bubble wall. Since the baryon number is generated at 
the wall in a manner sensitive to the wall velocity, these 
details are important for baryon number production. 

The tool used here is a hydrodynamical model where 
a scalar order parameter drives the transition [ p^ . 
The (non-perturbative) effective potential V{4)) for the 



scalar order parameter is parametrized by the latent heat 
L = L/Tc, the surface tension a = a jT^ and the corre- 
lation length Ic = IcTc- Non-equilibrium phenomena at 
the phase transition front are described by a simple fric- 
tion term, and the corresponding parameter 77 is deduced 
from a recent microscopic calculation Q| . 

The equations of motion follow from energy- 
momentum conservation and an equation for entropy 
production pO] . In an expanding universe with metric 



-dt^ + R{tf{dr^ + r^de^ + sin^eid^^) 



and in a spherically symmetric case, eqs. (5-7) of |Tl| ] 
become 
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Here e is the energy density and w = e + p the enthalpy 
density, and we have defined E = ej, Z = wj'^v. The 
velocity v above is the coordinate velocity dr/dt of the 
fluid, the physical velocity being Rv. The scale factor 
R changes very little during the period of interest, so 
the difference is unimportant, but below v will refer to 
the physical velocity. The main effect of the expansion 
is the cooling of the universe, where the small change is 
significant because we are very close to Tc. H = R/R is 
the expansion rate of the universe. 

For the thermodynamical properties of the phase tran- 
sition we use four sets of parameters, suggested by lattice 
simulations The main difference with respect to per- 
turbation theory is that the surface tension a gets small 
if the transition is weak. As the first set we take the 
Standard Model with an unrealistically small Higgs mass 
rriH — 51 GeV (case C): 



L = 0.124, a = 0.0023, Ic 



(4) 
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A second set is uih = 68 GeV in which case the transition 
is very weak but stiU of first order (case D): 



L = 0.08, a = 0.0002, L 



10. 



(5) 



Both of these cases lead to a transition too weak for 
baryogenesis. As a third case we consider a stronger tran- 
sition (case B): 



0.3, <T = 0.01, Ic 



(6) 



This transition corresponds to a tree-level Higgs mass 
parameter to^ ~ 50 GeV and is close to being strong 
enough for baryogenesis |^ (we do not consider the pos- 
sibilities proposed in |12]). In addition, the values in 
eq. (H) could be achieved in MSSM for a realistic pole 
Higgs mass mn ^ mw p3[ . A still stronger transition 
with ~ 45 GeV would give L ~ 0.5, a ~ 0.02 in 
MSSM. However, to display more clearly the parameter 
dependence, we take rather eq. (^ with just a larger sur- 
face tension a — 0.02 (case A). 

The dependence of the results on the parameters is 
roughly the following. The supercooling 1 — r„ (T„ = 
Tn/Tc) and the distances In/tn of the nucleated bubbles 
{tn = is the Hubble time, or Hubble length) are 

proportional to a^l'^ jh. The velocity of the bubble wall 
during the rapid growth is proportional to a^l"^ . Reheat- 
ing to Tc takes place if 1/ <^ Sir'^/'', and the velocity during 
the period of slow growth is proportional to (t'^/^/L^, pro- 
vided that this number is small enough. Hence a small 
tJ indicated by lattice simulations tends to make super- 
cooling and velocities small and reheating likely. 

The friction parameter r\ is determined from the mi- 
croscopic analysis in In general, the friction term is 
non-local and the coefficient r\ depends on jlij. How- 
ever, we do not here study the microscopic structure of 
the phase transition front and hence t] is only used in 
parametrizing the entropy production. Consequently, we 
may use a simpler 0-independent (but m//-dependent) t]. 
After determining L, ct, Zc and the nucleation temperature 
T„ corresponding to the effective potential used in |^ , we 
search for the value of r\ producing the same velocity as 
in Table 1 of @. We find that 77 - (0.04 - ^.\)T^ for the 
Higgs masses studied. For completeness we also inspect 
more general values of r\ below. 

For our four parameter sets A-D, we determine T„ and 
the distances Inj^H between the nucleated bubbles fl^ j 
with the following simple procedure, giving better accu- 
racy than the classical calculation. The parameters L, 17, 
If. fix uniquely a quartic potential V{(^\ and for this po- 
tential the nucleation action is known analytically with 
good accuracy [0. We then solve for T„ and In/tn nu- 
merically from eqs. (4.5), (4.8) in By comparing with 
nucleation calculations for more complicated functional 
forms of V{(j)), we have found that Ic should be inter- 
preted as the correlation length in the broken phase. 
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TABLE I. The properties of the transition. In case A 
there is no stage of slow growth, hence no foq, t^siow are given. 

After the nucleation the bubbles start to grow. We 
have studied the bubble growth with a spherically sym- 
metric hydrodynamical code, which solves eqs. ([l]"|l). We 
present results for 13 different sets of parameter values 
(cases A-D with different values of 77) , see Table I. 

The initial stages of bubble growth are not affected by 
the expansion of the universe or the existence of neigh- 
boring bubbles. The bubble wall accelerates and reaches 
a stationary velocity Wwaii- The velocity depends on the 
friction at the wall, parametrized by 77. If the super- 
cooling is relatively large, Uwaii covers most of the range 
from to 1 as i] is varied, and the configuration goes 
through the different modes described in ||l^. For small 
supercooling, only deflagrations are possible and there is 
a maximum wall velocity fmax p8[ . We plot fwaiK??) in 
Fig. Thus in case D the bubble will be a slowly grow- 
ing deflagration bubble. In cases B and C the bubble has 
to be a deflagration, but velocities closer to the speed of 
sound could be achieved. In case A detonations are pos- 
sible. Case D takes longer to compute, so we show Wwaii 
for 3 values of i] only. 

After the period of stationary growth, the bubbles col- 
lide with the compression waves from other bubbles. Our 
grid length corresponds to r„ = Zri/2, and we use reflec- 
tive boundary conditions at the (spherical) outer edge to 
represent collisions with the neighboring bubbles. This 
is an unrealistic geometry, since it makes the collision 
simultaneous at all points of the bubble wall, preserv- 
ing the spherical symmetry. Since the nucleation process 
places the bubbles randomly, in reality every bubble has 
a different collision geometry. The results presented here 
give a qualitative idea of the real events, although the 
aspherical features and transverse motions are missing. 

The bubble radius f = r/rn as a function of time for 
11 of the 13 cases is shown in Fig. |^. A fraction fv = 
of the volume is in the broken phase (bubble) and the 
transition is completed when f = 1. The microscopic 
(bubble nucleation radius) and the macroscopic (distance 
between bubbles) scales in the problem differ by some 9 
orders of magnitude. For practical reasons we have set 
the scales closer to each other on the computer. There- 
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FIG. 1. The wall velocity as a function of rj. In case A the 
wall is a detonation front for rj < 0.0077Tc. Otherwise the 
wall is a deflagration front. The black diamonds correspond 
to the 13 runs discussed in the text. 

fore the duration of the initial acceleration stage appears 
greatly magnified in the figure. 

We discuss the fastest- growing bubbles first. In case 
AO the bubble wall is a detonation front. No informa- 
tion is transmitted ahead of a detonation, and therefore 
the wall propagates undisturbed until it collides with the 
neighboring bubble wall. The phase transition is com- 
pleted and all the resulting fluid motion takes place in 
the broken phase. 

In all the other cases the bubble wall is a deflagration 
front, preceded by a compression wave in the symmet- 
ric phase. The bubble wall slows down as it meets the 
compression wave from the neighboring bubble. In most 
cases studied the compression wave is strong enough to 
stop the wall and push it back, compressing the bubble. 
As the compression waves move back and forth between 
the bubbles, the bubble volume oscillates. In a realistic 
geometry the bubble is deformed as different parts of the 
bubble wall move with different velocities. The typical 
amplitude of the motion should, however, be the same 
as obtained here. Part of the compression wave is trans- 
mitted through the wall into the broken phase, causing 
partial superheating. In some cases a region inside the 
bubble may be heated so much that the broken phase 
cannot exist at this temperature, and a droplet of the 
symmetric phase is temporarily formed. 

Further history of the bubble depends on the mag- 
nitude of the latent heat relative to the degree of su- 
percooling. If the latent heat is large enough to re- 
heat the universe back to the critical temperature, both 
phases can coexist at Tc- The initial volume fraction 
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FIG. 2. The bubble radius as a function of time. For clarity 
the cases A2 and B2 are omitted. 

of the broken phase at equilibrium (ignoring the ex- 
pansion of the universe since bubble nucleation) is then 
fy = f^^ = €c{l~T^)/L, where ec is the energy density of 
the symmetric phase at Tc- The bubble radius oscillates 
around fcq. The initial amplitude depends on Wwaii, as 
deflagrations close to the speed of sound are preceded by 
strong compression waves which cause large oscillations. 
A smaller f^q also leads to larger oscillations as the com- 
pression waves travel longer between the bubbles. The 
oscillation amplitude decreases with time and the bubble 
settles to a long period of slow growth. 

In case A the nucleation temperature is so low that 
the universe does not reheat to Tc. The transition is 
completed after at most one oscillation. The duration of 
the phase transition At is given essentially by the bub- 
ble separation and the wall velocity. At = r„/i;waii. In 
case Al the wall changes nature from a deflagration to a 
detonation after the collision. 

In cases B, C, and D, the universe reheats to Tc and 
the bubble radius oscillates around the equilibrium value 
fcq, given in Table I. 

The period of slow growth is paced by the expansion of 
the universe, which slowly makes space for the released 
latent heat. This takes a much longer time than the pre- 
ceding stages, see Fig. ||. To first order in L/ec, the dura- 
tion of the phase transition is At/tH = L{1 — fpq)/(4ec), 
assuming that in the symmetric phase p oc T^. The wall 
velocity Ugiow during the slow growth is geometry depen- 
dent, but for the spherical bubble geometry we can esti- 
mate it by fsiow ^ (1 — feq)(r„/At). The values of At/tn 
and fsiow are displayed in Table I. For case D Wsiow is 
much less than for the other cases, because the bubbles 
are much closer to each other. 
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FIG. 3. The bubble history showing the stages of fast and 
slow growth. The upper line is case B3 and the lower one CI. 



At some point during the slow growth the remaining 
regions in the symmetric phase puU themselves to spher- 
ical droplets, resulting in a geometry inverse to the one 
used above. We have done runs with an inverted ge- 
ometry, i.e., the shrinking symmetric phase at the center 
surrounded by the broken phase, starting at a static equi- 
librium situation at T ~ Tj,. The droplet evolves first to- 
wards the weak deflagration similarity solution described 
in and the end stages are as found in ITT] . 

The phenomena discussed might have physical signif- 
icance, e.g., for baryon number generation. Indeed, the 
velocity dependence of baryon number production may 
be such that the slowing down of Wwaii through colli- 
sions with compression waves, and the resulting oscil- 
lations, have a favourable effect on the baryon number 
produced It should be noted, however, that for 

stronger transitions (D^B) where the baryon number is 
not washed out afterwards, the fraction of space where 
the transition proceeds through slow growth is smaller. 

Finally, let us note that there remain some important 
open questions concerning the real-time history. In par- 
ticular, while the bubble wall appears to be stable during 
the period of fast growth ||2^ , the period of slow growth is 
different in character and the stability properties remain 
open. A related problem is the appearance of turbulence, 
which might affect, e.g., magnetic fields 
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